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This  technical  note  describes  a calculation  procedure  for  a pre- 
diction model  presently  in  use  at  the  Naval  Undersea  Warfare  Center, 

San  Diego  Division.  This  note  is  not  to  be  considered  as  an  official 
NUWC  report.  Its  purpose  is  to  document  the  existing  model  to  the  ex- 
tent required  by  Navy  project  offices  for  whom  predictions  have  been 
made . 
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his  note  describes  the  investigation  of  an  expression  to  be  used  in 
bottom  bounce  propagation  models  presently  in  use  at  NUWC  Code  D556. 
Recently  an  expression-1  was  developed  to  predict  the  acoustic  absorption 
coefficient  of  seawater  from  the  frequency,  temperature,  pressure  and 
salinity.  This  expression  was  developed  from  a combination  of  existing 
equations  in  order  to  provide  a single  prediction  expression  valid  under 
all  operating  conditions . 2 


In  the  past  many  methods  of  calculating  the  absorption  coefficient 
have  been  used.  These  methods  include  the  use  of  the  surface  temperature 

y . J 

and  50  F in  empirical  expressions.  This  paper  shows  that  it  is  necessary 
for  the  sake  of  accuracy  to  use  the  temperature  averaged  over  the  path  of 
the  propagating  ray  and  to  include  the  pressure  term. 

SUMMARY  AND  CONCLUSIONS: 

The  acoustic  absorption  coefficient  (a)  of  sea  water  for  any  particular 
point  is  a function  of  the  temperature  (T)  and  the  pressure  (P)  at  that 
point,  as  well  as  the  frequency  (f)  and  the  relaxation  frequency  (f^). 
Several  methods  of  averaging  a over  temperature-depth  profiles  were  compared 
for  accuracy  and  ease  of  computation.  The  difference  between  the  absorption 
calculated  using  depth-averaged  temperature  and  that  using  the  surface  tem- 
perature or  the  bottom  temperature  is  shown  in  figure  3.  The  temperature 
can  be  averaged  in  Procedure  RAYTRACE*  by  using  the  existing  values  of  depth 
and  time  increment  with  five  new  ETRAN  statements. 


^Procedure  RAYTRACE  is  a computer  program  which  computes  raypath  parameters 
using  Snell's  Law  and  a velocity  profile  comprised  of  constant  gradient  sections. 


The  absorption  loss  results  compare  favorably  with  those  calculated  using 


Errors  due  to  deletion  of  the  pressure  term  in  a calculations  can  be 


significant  at  large  depths  (>2x10  ft)  and  higher  frequencies  (>3kHz) 


Neglecting  the  increase  in  density, the  static  pressure  is  a linear  function 


of  depth  and  can  be  averaged  in  Procedure  RAYTRACE  with  the  addition  of  three 


new  statements 


PROCEDURE 


The  expression  for  the  acoustic  absorption  of  sea  water1  is 


02685f2 


1-6.  51+x10_4P 


dB/kyd 


32.768+f3  \l+32.768/f3 


where  f_  is  calculated  as  below 


( 6T+118 ) / ( T+21 3 ) 


The  average  temperature  is  calculated  from 


where  T is  the  temperature  during  the  time  period  t and  t is  the 


total  time  during  the  averaging  process.  The  average  relaxation  frequency 


same  manner 
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The  coefficient  of  absorption  was  calculated  separately  using  the 
average  temperature,  the  average  relaxation  frequency,  the  surface  tem- 
perature, and  the  bottom  temperature.  The  coefficient  was  also  determined 
for  each  break  in  the  temperature-depth  profile  and  averaged  to  provide  a 
physically  accurate  value  for  comparison. 

The  above  calculations  were  made  for  several  temperature-depth  profiles 
derived  from  bathythermograph  and  Nansen  cast  data  (figure  l)  assuming  that 
the  time  spent  at  each  temperature  is  proportional  to  the  depth  interval  of 
that  temperature.  Calculations  were  then  made  for  a general  profile  (figure  2) 
using  time  intervals  calculated  by  Procedure  RAYTRACE,  which  gives  the  actual 
time  spent  by  the  ray  in  each  depth  interval. 

The  values  of  a as  a function  of  frequency  resulting  from  the  applica- 
tion of  these  averaging  techniques  to  the  general  profile  are  shown  in  figure  3- 
It  is  seen  that  the  value  determined  through  temperature  averaging  is  nearly 
equal  to  the  corresponding  value  determined  by  averaging  the  coefficients  for 
each  temperature  along  the  profile. 

The  differences  between  the  absorption  loss  at  a range  of  50  kyd  de- 
termined by  the  surface  and  the  bottom  temperatures  are  shown  in  figure  U. 

Results  are  shown  for  several  operating  frequencies.  Figures  5,  6 and  7 show 
a comparison  of  a determined  using  the  surface  temperature,  the  bottom 
temperature,  and  the  average  temperature  for  several  varied  temperature 
profiles  at  a reference  distance  of  one  kiloyard. 

**  H 

The  variations  in  the  absorption  loss  caused  by  pressure  are  shown  in 
figure  8 for  several  frequencies  and  depths. 

3 


The  temperature  and  pressure  averaging  procedures  described  in  this 
note  are  being  implemented  in  a modified  Procedure  RAYTRACE4  which  will 


be  documented  in  the  near  future.  A direct  application  of  this  technique 
is  in  the  Bottom  Bounce  Propagation  Mode.  Adaptations  are  also  being 
considered  for  the  Convergence  Zone  Mode. 

A table  is  inserted  at  the  end  of  the  report  listing  the  absorption 
as  a function  of  frequency  and  temperature  for  several  pressures.  Zero 
pressure  is  also  listed  and  will  coincide  with  values  previously  presented.5 
The  depth  in  feet  required  to  produce  the  pressure  is  also  listed. 
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FIGURE  2 - GENERAL  PROFILE 
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Comparison  of_the  Absorption  Loss  calculated  from  the  average 

temperature  (T)  with  that  calculated  from  the  bottom  temperature 

(T  ) and  the  surface  temperature  (T  ). 
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Figure  8.  Comparison  of  a calculated  with  and  without  pressure  term. 


TABLES  - ABSORPTION  COEFFICIENT  vs  FREQUENCY  (kHz)  AND 
TEMPERATURE  (°F)  FOR  VARIOUS  PRESSURES  (DECIBARS) 
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